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Abstract 

Most wind farms consist of horizontal axis wind turbines (HAWTs) 
due to the high power coefficient (mechanical power output divided 
by the power of the free-stream air through the turbine cross-sectional 
area) of an isolated turbine. However when in close proximity to 
neighbouring turbines, HAWTs suffer from a reduced power coeffi- 
cient. In contrast, previous research on vertical axis wind turbines 
(VAWTs) suggests that closely-spaced VAWTs may experience only 
small decreases (or even increases) in an individual turbine's power 
coefficient when placed in close proximity to neighbours, thus yielding 
much higher power outputs for a given area of land. A potential flow 
model of inter- VAWT interactions is developed to investigate the effect 
of changes in VAWT spatial arrangement on the array performance 
coefficient, which compares the expected average power coefficient of 
turbines in an array to a spatially-isolated turbine. A geometric ar- 
rangement based on the configuration of shed vortices in the wake of 

*The technology described in this paper is protected under both US and international 
pending patents filed by the California Institute of Technology. 
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schooling fish is shown to significantly increase the array performance 
coefficient based upon an array of 16 x 16 wind turbines. Results sug- 
gest increases in power output of over one order of magnitude for a 
given area of land as compared to HAWTs. 

1 Introduction 

Technology for extracting energy from the wind has become a topic of in- 
creasing political, economic and social interest in the last fifty years [HI2]- 
Nearly all current wind farms use horizontal axis wind turbines (HAWTs) 
due to their high power coefficient, defined as the mechanical power out- 
put of the turbine divided by the power of the free-stream air through the 
cross-sectional area of the rotor normal to the free-stream direction. How- 
ever, a major constraint on wind farms is the inter-turbine spacing, since the 
performance of each HAWT decreases with decreased spacing [3H6]. When 
a HAWT extracts energy from the free stream, the speed of the flow be- 
hind the turbine is reduced. A well-known aerodynamic result is that the 
maximum amount of power a rotor can extract is proportional to the cube 
of the free-stream speed [I]. Horizontal axis wind turbines downstream of 
other HAWTs have a reduced incoming flow speed and thus are capable of 
extracting significantly less power than that of a spatially-isolated turbine. 
Theoretical analysis, computational simulations, and experimental results in- 
dicate that decreasing the downstream spacing of a row of multiple HAWTs 
from ten rotor-diameters to five rotor-diameters can decrease the power gen- 
erated by roughly 40% 017118]. Consequently, most modern HAWT wind 
farms space turbines five to ten rotor-diameters downstream and three or 
more rotor-diameters laterally [HE]. 

An alternative paradigm for wind energy extraction is found in verti- 
cal axis wind turbine (VAWT) designs. Studies have shown that a single 
spatially-isolated VAWT has a significantly lower power coefficient and costs 
approximately a third more than a HAWT with comparable power output [I]. 
Current literature regarding the performance of VAWT wind farms is scarce 
when compared to HAWTs, yet some sources indicate that the typical spac- 
ing constraints on HAWT wind farms do not apply [TUHT2] . Furthermore, 
a recent U.S. Patent suggests that aerodynamic interactions between a pair 
of counter-rotating VAWTs contribute to a higher power coefficient per tur- 
bine [13]. Previous work by Rajagopalan, Klimas, and Rickerl suggests an 
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increase of up to 4% in the average turbine power coefficient of an array of 
VAWTs due to steam-tube contraction (induced higher flow speeds around 
neighbouring VAWTs) [12J. In contrast, Schatzle, Klimas, and Spahr show 
no change or a mild reduction in power coefficient for a pair of closely-spaced 
VAWTs (TO] . However, neither of these studies considered the possibility of 
alternating the directions of rotation of neighbouring VAWTs. 

For this work, we consider a VAWT array configuration following the ar- 
rangement of shed vortices in the wake of schooling fish. These shed vortices 
form a reversed Karman vortex street in the wake of the fish [TJ] . Previous 
work by Weihs P3] uses a potential flow model to analyze the wake effects 
of the shed vortices on a school of fish. Motivated by the demonstrated ben- 
efits of the reversed Karman vortex street on the propulsion of the schooling 
fish, we apply the same configuration and similar modelling tools to analyze 
VAWT arrays. 

We investigate the inter-turbine spacing effects of VAWT arrays through 
a semi-empirical model. Models are developed, given in section [21 to de- 
scribe a single- VAWT and arrays of VAWTs using parameters from the field 
measurements. The field measurements of an installed VAWT are presented 
in section |3j The results of fitting the model parameters to the field mea- 
surements and evaluating effects of array configuration on power output are 
given in section @] and discussed in section [5j 



2 Modelling 

2.1 Potential Flow Model 

In the interest of examining the effects of spacing and orientation of VAWTs 
in an array, a potential flow model is first constructed for a single turbine. 
The model consists of uniform flow with a dipole and a point vortex su- 
perimposed to represent one VAWT. The point vortex accounts for the flow 
rotation caused by the turbine and the dipole represents flow blockage due 
to the presence of the turbine in the flow. The complex velocity potential, 
W, and the velocity vector field, u, are given by: 

W = U 00 z-t^-\og(z)+fiz- 1 (1) 

v v ' 

One VAWT 

u = U{dW/dz)t - %(dW/dz)j (2) 
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where Uoq refers to the free-stream speed, T is the strength of the point 
vortex, /i is the strength of the dipole, i and j are the unit vectors in the x 
and y directions, respectively, and z = x + iy (x and y are real numbers). 
To create the velocity potential corresponding to an array of K VAWTs, the 
model in equation [TJ is summed over the VAWT positions in the array, 



K 



W = U 00 z + J2 



k=0 



-i— log (z -z k )+/i(z- z k y l 

Z7T 



(3) 



kth VAWT 



where z k corresponds to the location of the kth VAWT. The velocity field is 
obtained using equation |2j To account for a velocity-deficit in each VAWT 
wake, the velocity-deficit curves given by Hau [I] are applied: 



u (z) 



(l ( 4 ) 



where u* is the resulting velocity vector after accounting for wake-effects 
and £ w (z) is a spatial function representing the velocity deficit. The function 
£ w (z) is a normal probability density function (PDF) in the lateral direction 

and a beta PDF in the downstream direction. The term ^1 — £ w (z)j adjusts 

the magnitude of the downstream velocity to reflect the velocity deficit in the 
wake. The wake extends approximately six rotor-diameters downstream and 
20° bilaterally from the free-stream wind direction. This velocity deficit is 
included for each VAWT in the array. Accounting for the wake in this manner 
implies that u* is no longer a solution to the potential flow equations. 

In order to evaluate the expected power produced by a turbine in an array, 
we define a parameter representative of the power output. By considering 
a simple drag-based VAWT, the tangential force on each turbine blade is 
~ u|, where Ut is the velocity tangent to the blade motion. Furthermore, we 
assume that angular velocity of the rotor is ~ Ut- Putting these together, 
the model for the power output of a single, spatially-isolated VAWT can be 
represented by: 

P lso = C lso ^- I (u*-sfds (5) 
27r Jo 

where P iso is the power output of the turbine, p is the density of air, C iso is 
a non-dimensional parameter specific to a VAWT design, Rc is the radius of 
the wind turbine, s is the unit vector tangent to the contour of integration, 
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and C is the contour of integration that is a circle with radius Rc centered on 
the VAWT in the direction of rotation. The non-dimensional parameter Ci SO 
is included to allow for the variety of wind turbine designs and accounts for 
effects such as blade shape, height of turbine, and other factors that affect 
wind turbine performance. The same expression for Pj SO (with a different 
C iso ) is obtained for a lift-based VAWT model as the lift is also ~ uf and it 
is still reasonable to assume the angular velocity of the rotor is ~ u t . 

Having defined Pi SO we introduce the array performance coefficient, Cap-, 
to evaluate the effects of closely spaced turbines in an array: 

C AP = %^ (6) 

where P ar ray is the average performance parameter for an array. It is implic- 
itly assumed that all turbines in an array are homogeneous (i.e. same Ci SO ). 
For Cap < 1, the average VAWT in the array performs worse than an equiv- 
alent, spatially-isolated VAWT. However, if Cap > 1, the model suggests 
that the array configuration provides beneficial, inter- VAWT, aerodynamic 
interactions and increases the average power output of the array. 

Although comparing Cap is a means for understanding the flow physics, 
as a means to compare the performance of an array configuration in terms 
of power output per unit land area occupied by the VAWT farm, a power 
density parameter, Cpd, is defined: 

Cpd = (^) I = K ■ Cap^- (7) 

V ^array J J \ ^iso J ^array 

where A array is the land area of the entire wind farm, Ai SO is the land area 
of a single turbine, taken to be arbitrarily D 2 , and K is the number of 
turbines in the array. This parameter is most useful in the practical design 
of a wind farm and determining the appropriate spacing for maximizing the 
power output per a given area. 

2.2 Biolnspiration 

A mathematical model similar to that described in the previous section is 
used by Weihs [15J to study the wake effects of fish schooling arrangements 
on fish propulsion. Weihs' potential flow model uses uniform flow and super- 
imposes point- vortices at the location of shed vortices from fish as they swim, 
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depicted in figure 1(a) The location of these shed vortices can be described 



in the complex plane by 
A 



m,n 

2+ 

m,n 



-a/2 + 2na + i(4mc + b) 
a/2 + 2na + i((4m + 2)c + b) 



(8) 
(9) 



A„ 
A 



m,n 

2- 
m,n 



a/2 + 2na + i(4mc - 6) 

-a/2 + 2na + z((4m + 2)c - b) 



(10) 
(11) 



where A]^ n and A^" n correspond to positive (anticlockwise rotating) vortices, 
A]^ n and A^ n correspond to negative (clockwise rotating) vortices, 2a is 
the downstream distance between two vortices of the same line, 2b is the 
lateral distance between the two lines of shed vortices from a given fish, 
2c is the distance between two adjacent fish, and m and n are integers. 
Fish in a school typically swim with b/a ~ 0.3, depending on species, and 
c/b > 2 [15j. Incorporating equation [9] and equation [TT1 into equation [3j the 
resulting velocity potential for an array of VAWTs following this particular 
arrangement is 



W = U 00 z 



^ E^l En=l [ log(* - A%) + \0g{z - <+J 

- log(z - A^ n ) - \og(z - A^- n ) 



2^m=l 2^n=l 



iV 



[z - A 



+ (z - A 



2+ \-l 
m,n) 



+(z- A 1 -)- 1 + (z-A 



2- x-li 
m,n/ 1 



(12) 



2.3 Numerical Simulation 

equation [12] is evaluated numerically for a VAWT array configuration where 
M = N = 8, corresponding to an 16 x 16 array of turbines (256 turbines 
total). The values of Uoo, T, and fi are obtained from the field measure- 
ments, as described in section 13.11 The model is interrogated for the range 
of parameters: ID < a < 8D, -1.5D < b < 1.5D, and c = 2D, where D 
corresponds to the turbine diameter. The range of a is chosen to extend from 
turbines immediately downstream of each other to where the turbines have 
sufficiently cleared the downstream wake of each other. The range of b is 
chosen to avoid overlapping turbines. 
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(a) Configuration of shed vortices in (b) Configuration for the VAWT wind 

the wake of schooling fish. Adapted farm based upon shed vortices in the 

from Weihs [15]. Coordinate system wake of schooling fish. Uoo is in the 

moves with the school, the fish are positive i direction, 
headed in the negative i direction. 

Figure 1: Diagrams comparing the model using shed vortices from schooling 
fish by Weihs [TJ)] with the arrangement of VAWTs used in the current study. 
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3 Field Measurements 



Field measurements of a Savonius-type VAWT located in Glendora, CA were 
collected in spring 2009 in order to generate input data for the potential flow 
model. The turbine was a prototype model with rotor radius of ~ 0.75 m 
and rotor height of ~ 2 m. In order to characterize the flow around the 
VAWT, two cup-type flow anemometers were used. A reference anemometer 
consisted of an adjustable- length pole with an anemometer (InSpeed Vortex, 
with a range of 1.4 - 55.6 m s _1 and ±0.05 m s -1 accuracy) mounted on the 
end and connected to a dedicated computer that recorded the wind speed 
every second. The pole was fixed to the ground using guy wires for untended 
operation. The data from the reference anemometer was averaged over the 
entire visit to produce the free-stream speed during the field visit. 

The second anemometer was mounted to a second pole that was manually 
positioned at 18 locations around the VAWT, as indicated in figure [2j In 
addition to measuring wind speed, the second anemometer contained a wind 
vane (InSpeed E-Vane) that measured the direction of the incoming wind 
(±1° accuracy). The wind vane was wired to a computer data acquisition 
system (Measurement Computing USB-1208FS) which recorded the voltage 
output of the wind vane at 100 Hz. This was then down-sampled through 
averaging the 100 measurements in one-second to obtain 1 Hz recordings 
that synchronized with the wind speed data. The anemometer recorded 
wind data at each station for approximately three minutes. For our field 
measurements, both instrument setups were extended to the mid-height of 
the VAWT and were adjusted according to the underlying terrain (« 6 m). 
At the end of the measurements, the mobile instrument setup was moved near 
the reference setup in order to obtain the free-stream direction, which was 
assumed constant during the one-hour measuring period. This is a reasonable 
assumption as local weather stations report only westerly wind during our 
measurement time period. Data from both anemometers was synchronized 
for post-processed based on the recorded time-stamps. 

3.1 Data Assimilation 

The model requires the calculation of three parameters to describe a single, 
spatially-isolated VAWT: Uqo, T, and \i. To obtain from the experimental 
measurements, the average free-stream speed from the reference anemometer 
is used. Using this measured free-stream velocity, the values of V and /x are 
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Figure 2: The measured flow field from the field visit to the VAWT along 
with the flow field calculated from equation [TJ using = 3.00 m s -1 , 
T = 7.41 m 2 s -1 , and /i = 0.00 m 3 s _1 . The solid vectors indicate the time- 
averaged experimental velocities. The model's least-squares-fit velocities are 
superimposed as the dotted vectors. Thin lines correspond to the external 
frame and rotor of the wind turbine. 

determined from a least-squares fit of a single-turbine model (one point vortex 
and one dipole) to the measured VAWT flow field. 



4 Results 

The measured flow field from the field visit to the VAWT along with the flow 
field calculated from equation [TJ using the model parameters described in 
section 13.11 are given in figure |5J The measured average free-stream speed, 
Uoo, is 3.00 m s _1 and appears to be well-described by a normal distribution 
having a standard deviation of 0.86 m s _1 . The fitting of the measured flow 
field to the single- VAWT model, as given in equation [IJ by least-squared 
difference yielded a unimodal minimum within the parameter space of V and 
fi. The measured results are well-represented by a range of values given by 
r=7.41±3.82 m 2 s" 1 and fi = 0.00±0.77 m 3 s -1 . Despite the range of U^, 
T and [A, in this letter we only present results using = 3.00 m s -1 , 
r = 7.41 m 2 s _1 , and y. = 0.00 m 3 s _1 . 

The calculated Cap for the fish-schooling VAWT-array models, given by 
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equation 



12.31 is shown in figure 3 (a 



and using the range of values of a and b described in section 
There is a clear region for which Cap is greater 



than unity, corresponding to close-spaced pairs of counter-rotating VAWTs. 
The maximum Cap in this parameter space suggests the potential for inter- 
turbine interactions that yield more power output than spatially-isolated tur- 
bines, figure 3(b) shows the normalized power density for the fish-schooling 
VAWT-array models. Comparing figures 3(a) and 3(b) shows that the opti- 
mal array configuration in terms of power extracted per unit land area does 
not coincide with the largest Cap values, which implies improving the aver- 
age performance of individual turbines might not be the desired goal of wind 
farm design. In figure [3J the dashed lines show the permissible parameter 
values for schooling fish. As seen in figure 3(a), this region is below unity 



in the Cap results. Array configurations of turbines with 4x4 and 32 x 32 
turbines were also explored and show results similar to the presented result. 



5 Discussion 

A qualitative comparison between the flow fields from the VAWT model and 
the collected field data, as shown in figure [2], demonstrates that the current 
single- VAWT model captures the dominant flow characteristics for a single 
turbine. This verifies the model's ability to reduce a complicated flow field 
to a three-parameter model. However from the discrepancy between velocity 
vectors of the single- VAWT model and collected field data in figure [2], it is ev- 
ident that the wake of a spatially-isolated turbine is not accurately captured. 
These discrepancies motivated the introduction of an empirically-based wake 
from HAWTs equation HI which is considered to be an appropriate means 
of accounting for an expected velocity deficit in the wake when calculating 
Cap- 

Measurements were limited to the area immediately surrounding the 
VAWT. Since the dominant flow effects of the VAWT are most pronounced in 
the closer locations, the vortex and dipole contributions to the velocity field 
drop off as r -1 and r -2 , respectively, where r is the distance from the vortex 
or dipole center. This explains our model's ability to capture the measured 
flow field despite the aforementioned limitation. Furthermore, measurements 
farther from the VAWT would likely reflect ambient fluctuations to a greater 
degree more than VAWT aerodynamics. 

In order to understand the practical implications of this work for wind 
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(a) Cap for the fish-schooling VAWT-array model. 
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(b) Cpd for the fish-schooling VAWT-array model, normal- 
ized by the maximum C pd ■ 

Figure 3: Calculated Cap and Cpd for the fish-schooling VAWT-array model, 
given by equation [TSJ are shown. Array-spacing parameters, a, b, and c 
as shown in figure 1(b), are normalized by the turbine diameter, D. For 
the above figures, c/D = 2. The black dashed line corresponds to typical 
parameters of schooling fish, b/a = 0.3 and c/b > 2. 
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Hypothetical VAWT Farms 

CE&P (Bio), 10 kW, 3 m dia. 
Helixwind S594 (Bio), 5 kW, 1.2 m dia. 
Mariahpower (Bio), 1.2 kW, 1.2 m dia. 
WePOWER Falcon (Bio), 1.2 kW, 1.78 m dia. 

Operational HAWT Farms 

Biglow Canyon Wind Farm, Oregon 
Blue Sky Green Field, Wisconsin 
Desert Sky Wind Farm, Texas 
Horse Hollow Wind Energy Center, Texas 
Smola Wind Farm, Norway 
Tararua Wind Farm, New Zealand 
Wild Horse Wind Facility, Washington 
Wildorado Wind Ranch, Texas 



Figure 4: Log plot showing the expected array power density, total array 
power output over array area, of VAWT and HAWT wind farms. Hypo- 
thetical VAWT wind farms using geometry based on fish schooling (a=1.2D, 
b=0AD, c=2D, Cap=0.61) are shown as the circles with solid lines. Geome- 
try chosen based on high Cpd- A sampling of operational HAWT wind farms 
are shown as squares with dashed lines. There is a significant increase of over 
one order of magnitude in the array power density for VAWTs compared to 
HAWTs. 
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farm design, a comparison of the array power density (total array power 
output over array area) between VAWTs and HAWTs is made, as given 
by Cpd. For this analysis, we assume that the parameters, F and fi, are 
the same for all VAWT designs. Additionally, we restrict ourselves to the 
parameter space representative of schooling fish, and chose the array config- 



uration that maximizes Cpp, following this constraint based on figure 3(b) 



figure H] demonstrates that significantly higher power densities are obtained 
with VAWT farms compared to operational HAWT farms, with increases 
of over one order of magnitude. Although close spacing of VAWTs poten- 
tially leads to higher power outputs of the individual turbines (regions of 
Cap > 1 in figure 3(a)), the major advantage of VAWT wind farms is their 
significantly smaller inter-turbine spacing, compared to HAWT wind farms, 
without a substantial decrease in performance as demonstrated by figures 



3(a) and 3(b) 



The proposed array configuration adapted from the geometry of fish 



schooling provides an interesting arrangement for study, figure 3(a) indi- 
cates decreases in the average turbine power output if turbines are geo- 
metrically arranged following the same configuration of shed vortices in the 
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wake of schooling fish. However the increase in array power density of close- 
spaced arrays overcomes these mild reductions in expected power output, as 
seen in figure HI If one allows for a broader range of parameter space than 
seen in schooling fish, there is the opportunity for significant improvements 
(Cap > 1) in the power output of these turbines. This is not entirely sur- 
prising due to differences in the criteria of beneficial arrangements. The fish 
aim to align themselves to optimize their forward propulsion whereas spatial 
configurations for turbine arrays aim to maximize energy extraction. 

While only one case is presented here, that of a 16 x 16 VAWT array, two 
other finite arrays were explored with dimensions of 4 x 4 and 32 x 32 turbines. 
These did not show any significant deviation from the case presented, and 
thus explains their omission. For all turbine array sizes explored, the wind 
farm configurations corresponding to the highest values of Cap, shown in 



figure 3(a), consist of pairs of closely-spaced, counter-rotating VAWTs with 
their adjoining axis perpendicular to the free stream. A streamline anal- 
ysis shows that stream-tube contraction between pairs of counter-rotating 
VAWTs, and the concomitant flow acceleration, is primarily responsible for 
improvements in expected power output. This is consistent with Rajagopalan 
et al. [12], who studied Darrieus-type VAWT arrays and showed a 4% in- 
crease in power coefficient due to the same mechanism. From our results 
it appears as though vortex interaction is the dominating mechanism for 
the stream-tube contraction, rather than blockage effects from neighbouring 
wind turbines. However, with further study, it is possible that this conclusion 
is dependent on VAWT design or array configuration, thus allowing VAWT 
designers to consider blockage effects in future design iterations. 



6 Conclusion 

In this work we present the experimental data used to create a potential 
flow model to describe the flow field around a spatially-isolated VAWT. This 
single- VAWT model consisted of three parameters: free-stream speed, vortex 
strength, and dipole strength. A model for an array of VAWTs, based on 
the single- VAWT model, is developed and an empirical wake is introduced to 
calculate the expected power output of the array. To evaluate the potential 
benefits of a particular array configuration, an array performance coefficient, 
Cap, is defined to compare the average power output of an array configu- 
ration to that of a spatially-isolated VAWT. Configurations inspired by the 
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arrangement of shed vortices in the wake of schooling fish are used for in- 
terrogating Cap values with varied geometric parameters. Results show the 
possibility of Cap values up to 1.4, for the 16 x 16 VAWT array, when not 
restricted to parameters corresponding to actual schools of fish. The leading 
mechanism for this increase is stream-tube contraction and concomitant flow 
acceleration between counter-rotating turbines in close proximity. However, 
values of Cap are approximately or less than unity for configurations follow- 
ing the schooling of fish regardless of array dimensions considered. These 
configurations significantly reduced the land use for VAWT wind farms, re- 
sulting in array power density increases of over one order of magnitude com- 
pared to operational HAWT wind farms. 

Although the current model does not explicitly allow for complicated flow 
behavior, such as vortex shedding, turbulence, and three-dimensional effects, 
the current vortex and dipole model yields a flow field sufficient for assessing 
the performance of VAWTs. Furthermore, the field data is only from a small 
segment of time (one hour for the particular visit reported) and does not fully 
describe the variety of conditions in which VAWTs can produce power. It 
is of interest to conduct further measurements on a broad range of full-scale 
VAWT designs and over different operating conditions in order to refine the 
single- VAWT model. Similar measurements and considerations are required 
to refine the model for arrays of full-scale VAWTs. Such refinements include 
investigation of the dependence of T and /i on VAWT design, wind speed, 
and surrounding environment (e.g. turbulence and the atmospheric boundary 
layer). Including these effects may potentially impact the results in terms of 
expected power output and array power density. 
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